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We determine the chiral phase transition line in (2+l)-flavor QCD for small values of the light 
quark chemical potential. We show that for small values of the chemical potential the curvature 
of the phase transition line can be deduced from an analysis of scaling properties of the chiral 
condensate and its susceptibilities. To do so we extend earlier studies of the magnetic equation 
of state in (2+l)-flavor QCD to finer lattice spacings, aT = 1/8. We use these universal scaling 
properties of the chiral order parameter to extract the curvature of the transition line at two values 
of the cut-off, aT =1/4 and 1/8. We find that cut-off effects are small for the curvature parameter 
and determine the transition line in the chiral limit to leading order in the light quark chemical 
potential. We obtain T c (fi q ) /T c (0) = 1 - 0.059(2)(4)(^ 9 /T) 2 + 0(ji\). 
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I. INTRODUCTION 

Extending lattice QCD calculations to non-zero baryon-chemical potential or, equivalently, to non-zero net baryon 
number density is known to be difficult in general. However, important information on the QCD phase diagram can be 
deduced for small values of the chemical potential by using well established numerical techniques such as rewcighting 
PJ , analytic continuation 0, Q or Taylor expansion 0, Q ■ We will concentrate here on the latter approach. 

Not only do Taylor expansions of thermodynamic observables provide information on QCD thermodynamics at small 
but non-zero chemical potential, the expansion coefficients themselves also are sensitive indicators for critical behavior 
in the vicinity of the chiral phase transition at vanishing chemical potential. As the chemical potential couples to the 
quark number current, which does not break chiral symmetry, it acts to leading order like a temperature variable. 
Derivatives with respect to chemical potentials lead to susceptibilities which exhibit critical behavior similar to that 
of thermal susceptibilities @ . We will show here that a calculation of the leading order Taylor expansion coefficient 
of the chiral order parameter, which defines a mixed susceptibility, allows to perform quantitative studies of the phase 
boundary between low and high temperature phases of QCD close to (i = 0. 

At non-zero values of the chemical potential a phase boundary in the temperature and chemical potential parameter 
space of QCD is well defined only in the heavy quark limit or for vanishing quark masses. In the former case the 
phase transition line corresponds to the first order deconfincmcnt transition in the pure gauge theory. At infinite 
values of the quark mass this transition is independent of the chemical potential and defines a straight line in the T-\i 
plane. For a large range of quark mass values the transition line is not unique. It characterizes a region of (rapid) 
crossover in thermodynamic quantities and a pseudo-critical temperature extracted from these observables may differ 
somewhat, depending on the observable that is used. In the chiral limit, however, the transition line is again well 
defined. For sufficiently large strange quark mass it defines a line of second order phase transitions in the universality 
class of three dimensional 0(4) symmetric spin models JtJ. 

Taylor expansions, analytic continuation as well as reweighting techniques have been used to locate the crossover line 
T c (ij) in the T-/J, plane for small values of the chemical potential Q. These calculations, which have been performed 
for different flavors and various values of the quark masses, suggest that the curvature is small, i.e. T c (/i)/T c (0) 
decreases only by a few percent at /i/T ~ 1. However, so far most lattice QCD calculations performed to determine 
the transition line have been performed on coarse lattice^, i. e. lattices with only four sites in the temporal direction 
for which the lattice spacing in units of the temperature thus equals aT =1/4. Better control over the extrapolation 
to the continuum limit and the quark mass dependence of the transition line clearly is needed. 

So far studies of the transition line concentrated on its dependence on a single chemical potential, taken to be 
cither identical for the light up and down ([id) quarks or of opposite sign. The former is the light quark chemical 
potential, /i q — {ji u + Hd)/2 and latter is the isospin chemical potential, (ii — — /id)/2 for which direct numerical 
calculations are possible. In order to make contact to the situation met in heavy ion collisions [ToL [Til ] one eventually 
wants to analyze the influence of non-zero charge (/xq) and strangeness (/is) chemical potentials on the curvature 
of the transition line, i.e., one should allow for non- vanishing up-, down- and strange-quark chemical potentials. 
For small values of the chemical potential this is possible in the framework we will outline here. At present we will, 
however, restrict our discussion to the case of vanishing strange quark and isospin (electric charge) chemical potentials. 

We concentrate on an analysis of the phase transition line in the chiral limit (m« = = 0, m s > 0) where its 
dependence on fx q is expected to be largest. We will present a calculation of the critical line for small values of the 
light quark masses in the scaling regime of the finite temperature chiral phase transition. This allows to use scaling 
relations to extract the curvature of the phase transition line in the chiral limit of QCD. The scaling relations naturally 
relate the curvature of the critical line as function of [i q to the magnitude of a mixed susceptibility. 

We will perform our numerical calculations for (2+l)-flavor QCD keeping the heavier strange quark mass close 
to its physical value and decreasing the two degenerate light quark masses towards the massless limit. On coarse 
lattice with temporal extent N T = 4 we will make use of a recently performed scaling analysis fl3j of the chiral order 
parameter performed with an improved staggered fermion action. This study showed that the chiral order parameter 
is well described by a universal scaling function characteristic for a three dimensional, 0(N) universality class. As we 
are using a staggered fermion discretization scheme for our scaling analysis we expect that the transition in the chiral 
limit, performed at non-zero lattice spacing, is controlled by the 0(2) rather than 0(4) universality class. We thus 
will analyze our numerical results in terms of 0(2) scaling functions. We will comment on the application of 0(4) 
scaling relations later on. 

This paper is organized as follows. In the next section we will extend the scaling analysis of the chiral order 
parameter to lattices with temporal extent N T = 8. This provides the basic parameters needed for a calculation of 



1 An attempt to determine the transition line closer to the continuum limit has been presented in Q . 
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the curvature of the chiral phase transition line which will be discussed in Section III. We give a discussion of our 
results and an outlook in Section IV. 



II. MAGNETIC EQUATION OF STATE 



In the vicinity of a critical point regular contributions to the partition functions become negligible in higher order 
derivatives and the singular behavior of response functions will generally be dominated by contributions arising from 
the singular part of the free energy densitjo 

f(T,mi,m s ,n q ,fi a ) = f a (T,mi,m a ,fi q ,fi a ) + f r (T,m h m s , n q , /i 8 ) . (1) 

In addition to the temperature T, light (mi) and strange (m s ) quark masses we also allow for a dependence of the 
free energy density on the quark chemical potentials. Close to the chiral phase transition temperature at vanishing 
chemical potential the singular part f s will give rise to universal scaling properties of response functions. This has 
been exploited to analyze basic universal features of the QCD phase diagram close to criticality fl2j . 

The singular part of the free energy density depends on the parameters of the QCD Lagrangian, e.g. the quark 
masses, and the external control parameters, temperature and chemical potentials, only through two relevant cou- 
plings. These scaling variables, t and h, control deviations from criticality, (t,h) = (0,0), along the two relevant 
directions, which in the case of QCD characterize fluctuations of the energy and chiral condensate, respectively. To 
leading order the scaling variable h depends only on parameters that break chiral symmetry in the light quark sector, 
while t depends on all other couplings. In particular, t will depend on the light quark chemical potential while h 
remains unaffected by these in leading order, 



_ 1 fT-T r . 



(2) 



to \ T c 
h = 1 1711 

where T c is the phase transition temperature in the chiral limit and to, ho are non- universal scale parameters (as is 

T c ). While the combination zo = h\/^ S /to is unique for a given theory, the values of to and ho will change under 
rescaling of the order parameter Note also that they depend on the definition of the parameter introduced to 
control symmetry breaking, i. e. the fact that we choose the strange quark mass to normalize the symmetry breaking 
light quark mass parameter. For the rest of this chapter we will not need to refer any further to the contribution of 
chemical potentials to the reduced temperature t. We will come back to it in the next chapter. 

The singular part of the free energy, f s , is a homogeneous function of its arguments. This can be used to rewrite 
it in terms of the scaling variable z = t/h 1 ^ 35 as 

f s (t, h) = h 1 + 1 / 5 f s ( Z , 1) EE h 1 + 1 / 5 f s (z) . (3) 

where (3, 5 are critical exponents of the three dimensional O(N) univer sality class HU, P = 0.349 and S = 4.780 
for three dimensional 0(2) models and /3 = 0.380 and 8 = 4.824 for (9(4), respectively. All parameters entering the 
definition of t and h, i.e. to, ho and T c may depend on the strange quark mass, but are otherwise unique in the 
continuum limit of (2+l)-flavor QCD. Just like the transition temperature T c , however, also to and ho are cut-off 
dependent and will need to be extrapolated to the continuum limit. 

The universal critical behavior of the order parameter, M ~ df/dmi, is controlled by a scaling function fa that 
arises from the singular part of the free energy density after taking a derivative with respect to the light quark mass, 

M(t,h) = hV'faiz). (4) 

The scaling function /g(-z) is well-known for the 0(2) and 0(4) universality classes through studies of three dimensional 
spin models This so-called magnetic equation of state, Eq. 2J has been analyzed recently for (2+lVflavor QCD 
using an improved staggered fermion formulation (p4-action) on lattices with temporal extent N T = 4 [13J and light 
quark masses as small as mi/m s = 1/80, which corresponds to a pion mass that is about half its physical value. It 
could be shown that the chiral order parameter can be mapped onto a universal O(N) scaling curve and the scale 



2 For systems belonging to the 3-dimensional 0(2) or 0(4) universality classes this does not hold for the thermal response function (specific 
heat) as the relevant critical exponent a is negative in these cases. 
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parameters to, ho, T c could be extracted. As the calculations had been performed with staggered fermions the scaling 
analysis has been performed by comparing results with the magnetic equation of state for a 0(2) universality class 
rather than 0(4) as one should find in the continuum limit for two massless quark flavor. However, as has been argued 
in [l3[ both scaling curves are similar in the limited range of z values, where this scaling analysis has been performed, 
and the scaling analysis could have been performed with 0(4) scaling functions as well. 

At least on these coarse N T = 4 lattices violations of scaling have been found to be small also for physical values 
of the light quark mass, i.e. mi/m s ~ 1/27. Of course, it is to be expected that the scale parameters, extracted 
on coarse lattices with temporal extent N T = 4, are subject to cut-off effects. We therefore extend the analysis of 
Ref. [IH to smaller lattice spacings. We perform calculations on lattices with temporal extent N T = 8. We follow 
here the discussion presented in Ref. [l3j and introduce two order parameters, that are multiplicatively renormalized 
by multiplying the chiral condensate with the strange quark mass, but differ in handling additive divergences, linear 
in the quark masa^l, 

M b = N*mS4>)i > 

M = N*m s - — {W)s) ■ (5) 



We use in our study data for the chiral condensate which have been collected by the hotQCD [lTj and RBC-Bielefeld 
[l8| collaborations in their studies of the (2+l)-flavor equation of state on lattices of size 32 3 x 8 as well as in the 
course of analyzing the transition temperature in (2+l)-flavor QCD [lj|. These calculations have been performed 
with the p4-action [2(J for three different quark mass ratios mi/m a = 0.2, 0.1 and 0.05, respectively. We use subsets 
of these data samples which cover a small temperature interval close to the transition region, but covering also the 
region of the transition temperature in the chiral limit. For the smallest quark mass ratio this includes a set of 16 
values of the gauge coupling which cover a temperature interval 0.95<T/T C <1.12 (3.48 < 6/g 2 < 3.545). Typically 
(20.000-30.000) trajectories of 0.5 time-units have been generated for each set of quark masses and gauge couplings. 
The temperature scale used in our scaling analysis is based on calculations of the scale parameter ro that has been 
determined from calculations of the heavy quark potential performed on lattices of size 32 4 [13, EH . 

The basic approach for the scaling analysis on N T = 8 lattices is identical to that described in Ref. [l3|. However, 
as the calculations on the N T = 8 lattices have not been performed for as small light quark masses as in the N T = 4 
analysis fl3j . where the smallest ratio was mi/m s = 1/80, we did not perform a separate analysis for the determination 
of scaling parameters in the small mass regime and a determination of scaling violating terms for larger values of the 
quark masses, as it has been done in Ref. [l3|. We perform a simultaneous analysis of data obtained for all three 
quark mass ratios and include scaling violating regular terms in the ansatz for scaling fits, 

M(t, h) = h 1/s f G (t/h 1/l3S ) + a t ATH + b x H + b a H 3 . (6) 

with AT = (T — T c )/T c and H = mi/m s . This includes all the scaling violating terms also used on N T = 4 lattices. 
In our final analysis, however, we will set & 3 = 0. 

We verified that this approach, applied to the N T = 4 data set, and restricted to the same mass range available 
now on N T = 8 lattices, i.e. 1/20 < mi/m s < 1/5, leads to results compatible with our earlier findings. We extract 
values for to, ho and T c which are similar to those obtained previously. In fact, results for fit parameters obtained 
from the analysis of different order parameters, M and turn out to be in even better agreement. Results of this 
new scaling analysis for the N T = 4 data set is shown in the upper half of Fig. [1] All fit parameter are summarized 
in Table |U 

Using the approach described above for the analysis of our N T — 8 data set we find good agreement of fit parameters 
extracted from an analysis of Mf, and M, respectively. Results of this scaling analysis are shown in the lower half 
of Fig. [TJ All fit parameter are summarized in Table HI As already noted in the analysis performed on lattices with 
temporal extent N T = 4 jl3| . we observe also for N T = 8 that scaling violations are small for mi/m s < 1/10. This 
confirms that physical quark mass values, corresponding to mi/m s ~ 1/27, are in the scaling region. 

The constants ho, to determined for N T = 8 take on values different from those for N T = 4. The invariant 
combination of scale parameters, 

zo = h 1 o /0S /to = zo(m s ) + O(a 2 ) , (7) 

changes by about 50% which shows that its continuum extrapolation is not yet possible. This also is the case for to 
and ho separately. 



3 At finite values of the cut-off these terms are. of course, finite and may be viewed as a specific contribution to the regular part that 
will not alter the scaling properties for sufficiently small values of the quark mass. 
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FIG. 1. Fit of the 0(2) scaling function to numerical results for the subtracted order parameter M (right) and the non- 
subtracted light quark condensate Mb (left), both for N T = 4 (top) and N T = 8 (bottom). The fits include an ansatz for 
violations of scaling as discussed in the text. Shown are results for rai/m s < 1/5. 



When comparing results obtained for N T — 4 and N T = 8 one also has to take into account the dependence of the 
scale parameters on the strange quark mass. In fact, as the scaling analysis has been performed at bare strange quark 
mass values m s , fixed in lattice units, the corresponding physical value in the chiral limit at t = is only determined 
a posteriori, once T c has been determined. It turns out that the physical values of the strange quark mass in the 
N T = 4 and 8 calculations differ at T c by about 10%. One may account for this mismatch by rcweighting the results 
for the chiral condensates in the light and strange quark masses [23j | . However, we will not attempt to do this here. 

The main outcome of the N T = 8 scaling analysis, aside from confirming the good scaling properties of the order 
parameter at a twice smaller value of the lattice spacing, is a determination of the scale parameters and the transition 
temperature T c in the chiral limit, needed in the definition of the scaling variable z, i.e. the determination of io, ho 
and T c . We summarize these results in Table HI In the next section we will make use of these scale parameters to 
determine the curvature of the phase transition line for small values of the quark chemical potential. 



III. CURVATURE OF THE CRITICAL LINE 

As outlined in the beginning of the previous section at leading order the light quark chemical potential only enters 
the reduced temperature t, as introduced in Eq. [2J Also at non-vanishing values of the quark chemical potential 
the phase transition point is located at t = 0. The variation of the transition temperature with chemical potential 
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0.0037(2) 


0.0022(3) 


194.5(4) 






6.8(5) 




M 


0.0048(5) 


0.0048(2) 


195.6(4) 
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M b 


0.00407(9) 


0.00295(22) 


194.9(2) 


3.8(21) 


2.1(1) 


7.5(3) 




M 


0.00401(9) 


0.00271(20) 


194.8(2) 


11.2(21) 


-2.4(1) 


7.2(3) 
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M b 


0.00271(21) 


0.00048(9) 


174.1(8) 


-10.1(16) 


3.3(5) 


3.8(5) 




M 


0.00302(22) 


0.00059(10) 


175.1(8) 


-0.9(15) 


-4.6(4) 


3.8(4) 



TABLE I. Scale parameters determined from the scaling fits on lattices of temporal extent N T — 4 and 8. In column 6 and 7 
we list the couplings for the leading scaling violating corrections. The last column gives zo = h\ /to. We give the results for 
parameters entering the definition of scaling functions for Mb and the subtracted order parameter M as defined in Eq. [S] Only 
the former has been used in the analysis of the mixed susceptibilities. Note that fits including regular terms, give consistent 
determinations of the parameters of the scaling functions determined from Mb and M, respectively. 



therefore is parametrized in terms of the constant n q introduced in Eq. [2] 

T4^ = i-«,(%)\o^). (8, 

To determine the chiral phase transition line in the T-fi plane we thus need to determine the proportionality constant 
K q . This is, in fact, the only left over free parameter in universal scaling functions that needs to be determined. 
All other parameters (to, ho, T c = T c (/i q = 0)) have already been determined in the scaling analysis of the order 
parameter discussed in the previous section. 

The constant n q can be determined by analyzing the dependence of the chiral condensate on the light quark 
chemical potential. Of course, at vanishing light quark mass one would simply determine the temperature at which 
(iptp)i vanishes. At non-zero but small values of the quark mass this information is encoded in scaling functions. To 
extract information about the dependence of the scaling variable t on n q it suffices to consider the leading order Taylor 
expansion coefficient of the chiral condensate, 



(9) 



where 



Xm, g _ d 2 W)z/T 3 _ d Xq /T* 
T d{^ q /Tf dmi/T ■ 1 ' 

The mixed susceptibility Xm,q is proportional to the leading order coefficient of the Taylor expansion of the chiral 
condensate, which has been introduced in [2lL l22j . It may also be viewed as the quark mass derivative of the light 
quark number susceptibility (Xq)- Details of its definition in terms of inverses of the staggered fermion matrix and its 
derivatives with respect to the quark chemical potential are given in Appendix A of Ref . [22} • We have summarized 
the formulas relevant for our current analysis in an Appendix. 

In the masslcss limit the chiral order parameter vanishes at T c and varies as M ~ (— t)@ . Its derivative with respect 
to t thus will diverge at T c like t p ~ x . The same singular behavior will thus show up in a derivative of the chiral 
condensate with respect to temperature as well as the second derivative with respect to fi q /T. The pre-factors of the 
singularity in dM/dT and d 2 M/d(/i g /T) 2 , however, will differ by a factor 2n q T c . We will make use of this relation to 
determine the curvature of the critical line at \i q = 0. 

In the vicinity of the critical point the mixed susceptibility can be expressed in terms of the scaling function 

f' G {z) EE df G ( Z )/dZ, 



T torn 



The sca ling function f' G {z) is easily obtained from fa(z) by using the implicit parametrization for the latter given 
in Ref. [16|. We also note that Xm,q diverges as function of the light quark mass at t = 0, i.e. at the chiral phase 
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N r 


mi/m s 




2 /J £ 

X /dot 


4 


1/10, 1/20 


0.0598(26) 


3.5 




1/40, 1/80 


0.0573(29) 


1.5 


8 


1/20 


0.0559(35) 


0.4 


4, 8 


all 


0.0591(17) 


2.1 



TABLE II. Determination of the curvature of the critical surface of the chiral phase transition in (2 + l)-flavor QCD as function 
of the light quark chemical potential \i q . The table summarizes fits performed separately for two lighter and two heavier quark 
mass sets as well as the combined data set. 



transition temperature. In contrast to the chiral susceptibility, \ m ~ dM/dmi, which stays finite in the chiral limit 
only for t > 0, the mixed susceptibility is finite for all t ^ 0; for t < it scales like (— t) 13 ^ 1 while for t > it behaves 
like /it" 1 "'*'. 

For small values of the light quark mass numerical results for the mixed susceptibilities Xm,q may be compared to 
the right hand side of Eq. lTTI Here all parameters that enter /q(z) are known and the only undetermined parameter is 
Kg. As we did for the analysis of the magnetic equation of state, we should consider the influence of scaling violations 
induced by non-zero values of the light quark masses on the determination of the curvature of the phase transition 
line. The leading quark mass corrections identified in the analysis of the magnetic equation of state will not contribute 
to Xm.q as they do not depend on the chemical potential. The first scaling violating term would arise from a regular 
term that gives corrections to the order parameter of the form M ~ a q H(fj, q /T) 2 . This would give rise to corrections 
to the scaling relation given in Eq. QT] 

t^-ey^x^ = 2Kq f G ( z ) + ^hil H i+a-P)/^ (12) 

T T z 

As we do not know the pre-factor a q we need to check in the analysis of the mixed susceptibilities whether corrections 
to scaling play a role. We note, however, that in the case of the magnetic equation of state the dominant corrections 
arise from the term b\H , which is ultraviolet divergent in the continuum limit. This term drops out in the analysis 
of Xm.q- Moreover, scaling violating terms are further suppressed by a factor ij( 1 ~' 3 )/' 315 = ( or ff0.34 - m 0(4) 

symmetric models) as the dominant scaling term itself is divergent at (t, h) = (0, 0). We thus expect scaling violations 
to be small. 

Using a subset of the data samples described in the previous section, we calculated the mixed susceptibility Xm.q 
on lattices with temporal extent N T = 4 for several values of the quark mass. For this analysis we used data sets 
separated by 50 trajectories. For the lightest quark mass ratio, mi/m s = 1/80, we selected 4 and for the three heavier 
quark mass ratios, mi/m s = 1/10, 1/20, 1/40, we choose 6 values of the gauge coupling in a narrow temperature 
interval close to the chiral phase transition temperature T c , i.e. —0.02 < (T — T c )/T c < 0.06. Typically this involved 
about 500 to 950 gauge field configurations per parameter set, except for the lightest quark mass ratio where we 
analyzed about 350 gauge field configuration. On each gauge field configuration we calculated the various operators 
necessary to construct Xm,q- Explicit expression for them are given in the Appendix (see also Appendix A of (22j|). 

The calculation of the various operators required inversions of the staggered fermion matrix with a large set of 
random noise vectors. We used 500 noise vectors on each gauge field configuration and constructed unbiased estimators 
for the various traces that need to be calculated. All these calculations could be performed very efficiently on a GPU 
cluster. 

Results obtained for the mixed light quark number susceptibility, Xm,q, on lattices with temporal extent N T = 4 
are shown in Fig. [2j We clearly see that Xm,q increases in the transition region with decreasing values of mi/m s . 

Using the scaling relation given in Eq. [TT]wc can re-scale the data and obtain a unique scaling curve. This scaling 
curve can be mapped onto the 0(2) scaling function /q(z) with a simple multiplicative rescaling factor, 2n q . The 
resulting scaling plot is shown in Fig. [3J To check for possible contributions from scaling violating terms we have 
analyzed the data separately for quark mass ratios mi/m s = 1/10, 1/20 and m;/m s = 1/40, 1/80. These fits agree 
within statistical errors. We then determine the curvatures n q from fits to the complete data set. Results of these fits 
are summarized in Table [TTJ 

The scaling analysis performed for the mixed susceptibility on lattices with temporal extent N T = 4 suggests that 
the determination of the curvature parameter k m can be reliably performed with quark masses m;/m s <l/10. This is 
in accordance with the scaling analysis of the order parameter itself, which we have discussed in the previous section. 
It thus seems to be safe to extract the curvature parameter also at smaller values of the lattice spacing, i.e. from 
our N T — 8 data set, by using the smallest quark mass ratio available there, mi/m s = 1/20. We have performed 
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FIG. 2. The mixed light quark number susceptibility as a function of the reduced temperature, (T — T c )/T c . Shown are results 
obtained at two values of the cut-off, N T = 4 (open symbols) and N T = 8 (filled symbols), and for several values of the light to 
strange quark mass ratio. 

calculations at five values of the temperature using gauge field configurations on 32 3 x 8 lattices generated by the 
hotQCD collaboration [lj|. For these parameter sets we have analyzed 300 to 600 gauge field configurations, which 
were separated by 100 trajectories. Again we used 500 noise vectors for the calculation of all relevant operators on 
each of the gauge field configurations. The result of this analysis is shown in Fig. [2] and Fig. [3] with filled symbols. As 
can be seen they agree well with results obtained on coarser lattices. 

When rescaling data obtained for Xm, q to the 0(2) scaling curve f' G {z) we need to take into account errors on the 
scaling parameters to and z (or ho). This leads to a 10% error for the determination of the curvature terms. 

Performing a combined fit to all results obtained for different quark mass values and lattice spacings we obtain 

Kq = 0.059(2)(4) . (13) 

This result for the curvature of the critical line is about a factor two larger than the reweighting results obtained 
m (2+l)-flavor QCD Q. It is however consistent with results obtained in calculations with imaginary chemical 
potentials. In fact it lies inbetween the 2-flavor Q and 3-flavor j24[ simulations performed with the standard staggered 
fermion formulation and also is consistent with results reported from (2+lVnavor simulations with imaginary chemical 
potential performed with the action used also in this study (p4-action) [25|. 

IV. CONCLUSIONS 

With this analysis we have established a systematic way to determine the curvature of the QCD phase transition 
line in the chiral limit for small values of the light quark chemical potential. We have determined the curvature for two 
values of the cutoff using an improved staggered fermion action (p4-action) . Within our present statistical accuracy 
we did not observe any significant quark mass dependence of the scaled mixed susceptibilities. The result observed for 
the curvature of the second order phase transition line in the chiral limit thus also is a good estimate for the crossover 
line at physical values of the light quark masses. 

Although the final result for the curvature term seems to show little cutoff dependence, one has to be cautious as 
the other three scale parameters that enter this analysis (to, hg and T c ) all vary significantly as the lattice spacing is 
reduced by a factor of two. Clearly more work is needed to extrapolate safely to the continuum limit We will do so 
in the future by repeating the analysis with a discretization scheme that also suppresses cut-off effects arising from 
taste symmetry violation in the staggered fermion action more efficiently (hisq action) . 

In our current analysis we have kept the strange quark and isospin chemical potentials equal to zero. A direct 
comparison to the situation met in heavy ion collisions on the freeze-out curve thus should be done with caution. 
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FIG. 3. The scaled mixed susceptibility as function of the scaling variable z = t/h 1 ^ 136 . The data are compared to the 0(2) 
scaling curve. The band shows a 10% error band on this curve which arises from statistical errors on the calculated observables 
as well as from the errors on the scaling parameters to and zo given in Table [U 



However, at least for large energies, i.e., small values of the baryon chemical potential experimental results for the 
freeze-out curve correspond to electric charge chemical potentials, which are more than an order of magnitude smaller 
than fiB- As susceptibilities obtained by derivatives with respect to strange quark chemical potentials rather than light 
quark chemical potentials are generally smaller, one may also expect that the curvature in the [is direction will turn 
out to be smaller. It thus seems that the curvature of the critical surface along the /i q ~ /is /3 direction is most relevant 
for a comparison of lattice QCD results with the experimentally determined freeze-out curve. The phenomcnological 
parametrization of the freeze-out curve given in [ll| yields Tf reeze (p B ) /Tf reeze (0) ~ 1 — 0.21 (2) (// ? /T) 2 + 0(^ q ). 
The curvature of the freeze-out curve thus is about a factor 4 larger than that determined here for the chiral phase 
transition curve. This suggests that the freeze-out curve may not follow the chiral phase transition or crossover line 
at non-zero values of the chemical potential. With increasing fJ- q /T the hadronic freeze-out seems to happen further 
away from criticality. At the largest value of the light quark chemical potential currently explored in the low energy 
scan at RHIC [26[, /i q /T ~ 1, the freeze-out temperature may be about 15% below the crossover temperature. 

Nonetheless, as pointed out above, one still needs to improve the current lattice calculations. Results closer to the 
continuum limit with further improved fermion discretization schemes are needed and one should also get control 
over the influence of non-vanishing strange quark chemical potentials in order to firmly establish the separation of 
the freeze-out curve from the chiral transition line as advocated above. 
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Appendix A: The mixed susceptibility \r, 



We summarize here the operators entering a calculation of the mixed susceptibility Xm,q introduced in Eq. llQI This 
susceptibility is proportional to the second order Taylor ex pan sion coefficient of the chiral condensate in terms of the 
light quark chemical potentials. Using Appendix A of Ref. |22j ) we introduce the expectation value of the light quark 
chiral condensate as 



T 3 



(Co) 



(Al) 



and obtain for the mixed susceptibility, 



T 



1 

m 



1 d 2 {^)jT 3 



nb 



dp 2 



((C 2 ) + 2 (CiTh) + (C Q T> 2 ) + (C V 2 ) - (Co) ((V 2 ) + (V 2 ))) 



(A2) 



where we have introduced the shorthand notation p = \x q a for the chemical potential expressed in units of the lattice 
spacing. Here C n and T> n denote n-th derivatives of the trace of the inverse fermion matrix (D) and logarithms of its 
determinant, respectively, 



r 



1 d^rD- 1 
I" 



dp" 



V n = 



1 <9"lndetL> 



dp n 



(A3) 



where the derivatives are defined with respect to the flavor chemical potential p = pf for the quark flavor /. The 
factors 1/4 and 1/2 arise because we define the mixed susceptibility as a derivative of the 1-flavor light quark chiral 
condensate with respect to the light quark chemical potential fi q which is identical for the two light flavor components 
of the fermion action, i.e. a flavor factor n/ = 2 arises only in derivatives of the logarithm of the fermion matrix. 

thus arc 



The operators needed to calculate x 

Co 
Ci 

Co 

V 1 
V 2 



m.q 
-X") 



D -^° D -^° D -i 



1 ( t\ i dD 
2 tr ( D ^p 
1 / / 1 d 2 D 
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